Density functional theory for Baxter's sticky hard spheres in confinement 
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It has recently been shown that a free energy for Baxter's sticky hard sphere fluid is uniquely 
defined within the framework of fundamental measure theory (FMT) for the inhomogeneous hard 
sphere fluid, provided that it obeys scaled-particle theory and the Percus-Yevick (PY) result for 
the direct correlation function [Hansen-Goos and Wettlaufer, J. Chem. Phys. 134, 014506 (2011)]. 
Here, combining weighted densities from common versions of FMT with a new vectorial weighted 
density, we derive a regularization of the divergences of the associated strongly confined limit. 
Moreover, the simple free energy that emerges is exact in the zero-dimensional limit, leaves the 
underlying equation of state unaffected, and yields a direct correlation function distinct from the 
PY expression. Comparison with simulation data for both the bulk pair correlation function and 
the density profiles in confinement shows that the new theory is significantly more accurate than the 
PY-based results. Finally, the resulting free energy is applicable to a glass of adhesive hard spheres. 

PACS numbers: 61.20.Gy, 82.70.Dd, 64.60.ah 



Colloidal suspensions arc known to exhibit a variety 
of collective phenomena such as crystallization, gel for- 
mation, glass transition and percolation [iHll- Unlike 
in molecular systems these properties can be measured 
not only in terms of bulk sampling such as x-ray diffrac- 
tion, heat capacity measurements, or mechanical prop- 
erties, but also on a microscopic level thanks to opti- 
cal single-particle observation techniques such as confo- 
cal microscopy. An ideal theoretical tool to study the 
structural properties down to the particle size is classical 
density functional theory (DFT) an excellent exam- 
ple application being the purely entropic hard-sphere sys- 
tem. Indeed, within the framework of fundamental mea- 
sure theory (FMT) [l-Q, for arbitrary external poten- 
tials both the inhomogeneous liquid and the crystalline 
phases are described by a simple and accurate free energy. 
However, a more realistic model of colloidal suspensions 
must account for attractive interactions. Such a model 
is provided by Baxter's sticky hard spheres (SHS) which 
interact via an infinitely narrow and deep attractive well, 
the limit being taken such that the second virial coeffi- 
cient remains finite @ . The model is particularly appeal- 
ing because it is soluble analytically within the Percus 
Yevick (PY) approximation. The PY results obtained 
by Baxter provide a framework to interpret experimen- 
tal phase diagrams of colloidal suspensions, including the 
liquid-liquid transition, crystallization, percolation and, 
to some extent, the glass transition, revealing the exis- 
tence of attractive and repulsive glasses Q . Quantitative 
interpretation of single-particle observations in colloidal 
suspensions rely on a theoretical approach capable of re- 
solving the microscopic structure of Baxter's SHS as ac- 
curately as does FMT for hard spheres. Previous DFT 
approaches relied on Taylor expansions of the free energy 
around the homogeneous bulk fiuid which use the PY re- 



sult for the direct correlation function while crudely ap- 
proximating higher order terms. Recently, using a set of 
weighted densities from FMT for the hard-sphere fluid 
[lo| . a first non-perturbative DFT was derived. When 
applied to density profiles near a planar wall, or the SHS 
confined between two parallel planar hard walls, the the- 
ory was shown to be a significant improvement over pre- 
vious approaches. However, based on experience with 
Rosenfeld's original FMT for the hard-sphere fluid, it is 
expected that the FMT for the SHS would exhibit an un- 
physical divergence for strongly conflned fluids, the most 
extreme test being the zero-dimensional (OD) limit. 

In this Letter, we recast the original FMT for the 
SHS, which is based on Kierlik's and Rosinberg's scalar 
weighted densities 0, 0, by combining Rosenfeld's vec- 
torial weighted density [5| with a new vectorial weighted 
density. In consequence, we identify the origin of the di- 
vergence and regularize it in the spirit of Schmidt et al.^s 
approach to Rosenfeld's FMT for the hard-sphere system 
Rather than deriving a cumbersome expression by 
conserving the underlying PY correlation function, we 
seek a simple density functional. When the stickiness 
of the particles is strong, the modified theory (exact in 
the OD limit) provides a better description of the pair- 
correlation function than does the PY-based FMT, while 
the accurate results at moderate stickiness are left vir- 
tually unchanged. For the SHS confined in a spherical 
cavity, the new theory is significantly better than its pre- 
decessor. Our free energy can be used to study a wide 
range of the structural properties of inhomogeneous col- 
loidal suspensions with short-ranged attractions. In par- 
ticular, this is the first free energy applicable to the glassy 
state where confinement due to caging is ubiquitous. In 
fact, it may overcome previous theoretical pathologies as- 
sociated with Baxter's adhesive hard-sphere treatment of 



2 



the glass transition. These problems were based on using 
the P Y direct correlation function alone [l2| , rather than 
the complete free energy provided here. 

Consider a fluid of monodisperse hard spheres with ra- 
dius R. We define the locally averaged packing fraction 
n3(r) ~ J dr' p{r')<d{R — |r — r'|) with weighted densi- 
ties (r) = Wi^^^^^ (r) — — Vn3(r). Rosenfeld's 
excess free energy density (measured in units of ksT) 
is then $hs = $i + $2 + 'I's with $1 = -712 ln(l - 
n3)/(47ri?2), $2 - [nj - 1x2 ■ n2)/[47ri?(l - rig)], and 
$3 = {nl - 3n2n2 • n2)/[247r(l - n^f] % While $hs 
has been extremely successful in describing an uncon- 
fined fluid, with an inhomogeneity caused, say, by a pla- 
nar wall, it is problematic when the fluid is highly con- 
fined (see and refs. therein). For the most extreme 
confinement, the density distribution becomes simply a 
delta-peak, i.e., p^^{r) = 7]S{r), where r = |r|, and rj 
is the packing fraction (0 < rj < 1). The weighted 
densities for this distribution are n.3^(r) — ry6(_R — r), 
"2^(^) = ri5{R — r), and n2^(r) = r]S{R — r)r, where 
r = r/r. The integration of $1 can be readily performed 
yielding J dr<I>i = 77 -|- (1 — 77) ln(l — 77) which is the ex- 
act result for a OD cavity holding an average number of 
particles 77 [llj . Examining $hs further for such a cavity 
shows that $2 vanishes (as it should) while $3 exhibits 
a strong negative divergence. This is due to the factor 
3 in front of the vector term in $3 which would have 
to be unity for its scalar and vectorial terms to cancel. 
However, the resulting expression for the excess free en- 
ergy does not yield the PY direct correlation function, 
thereby destroying the excellent experimental agreement 
of Rosenfeld's $hs for the HS fluid. Several routes have 
been devised to regularize $3. A particularly appealing 
approach consists of introducing a new weighted density 
77-2 = 7i2 — (n2 • n2)/7i2 that reduces to 7i2 in the bulk 
fluid (where the vectorial weighted densities vanish as 773 
is spatially uniform), and vanishes in the OD limit (where 
n2 ■ n2 = 772). Using 77.2, a regularized excess free energy 



density proposed by Schmidt et al. |11| is 
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772772 
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(1) 

From the properties of 772 it is immediately clear that 
$HS is exact for a OD cavity while the free energy of 
the bulk fluid is unaffected. However, the third term 
differs by order [(n2 • n2)/7i2]^ from that of $hSi which 
is very small for moderately inhomogeneous fluids. In 
fact, it affects only the direct correlation functions of 
order 4 and higher whereas the PY direct correlation 
function (which is second order) is preserved. Both un- 
confined and confined fluids are extremely well described 
by $Hs; accurately predicting solid-liquid coexistence in 
hard sphere systems whereas ^hs dramatically oversta- 
bilizes the crystal due to the negative divergence for 
strongly peaked density distributions. However, Tara- 



zona's modification of $3, which uses a tensorial weighted 
density [l3| . performs better with regard to the width of 
the density peaks in the crystal. 

How does the picture change as we introduce Baxter's 
attractive surface interaction into the system? Rosen- 
feld's set of weighted densities (i.e. 773, 712, and n2) is 
known to be insufficient for constructing any reason- 
able free energy density for the SHS fluid because they 
cannot yield the (5-function within the direct correlation 
function [lo| . However, when 773 and 772 are supple- 
mented by the weighted density 771 (r) = g^^'T-2(r)j 
the set is sufficient to construct a free energy density 
•i'sHS = ^us + *&s for Baxter's SHS fluid where $s = 
ni(/)i/i?-|- 772^2/(27ri?^) is uniquely defined by (a) requir- 
ing the density functional to yield the PY direct corre- 
lation function and (b) imposing consistency with scaled 
particle theory [f3|. The dimensionless coefficients 0i 
and (ho are functions of x — .^"^ , which in the bulk 

fluid is related to the packing fraction 7/ via x — 
They are obtained by solving (j)[{x) = —2ycr{x)/x and 
[x^(j)2{x)]' = jjcix)^ /2 — xya{x) + xy'^{x), where the inte- 
gration constants are chosen such that and 02 vanish 
for a; — >■ 0. Here j/o- = riy^/r, where y^ is the cavity func- 
tion at contact, and is obtained as the smaller of the two 
solutions to yl - {Ax + l2T)yc + 2a::(2 + x) = 0. 

In order to regularize the divergences in the OD limit 
we rewrite $3 using vectorial weighted densities so that 
the resulting functional i>s vanishes. This guarantees 
that $Hs -I- is exact for a cavity that can hold only 
one particle at a time, and thus the corresponding free 
energy is unaffected by the sticky interaction. We use 
the identity 77i(r) = 772(r)/(47ri?) + (V • V)773(r)/(87r) 
in $s, and for reasons that will become clear below, 
we split the term ni(f)i/R in half and substitute for 771 
in only one of the terms. Moreover, because $s is al- 
ways integrated on R^, we can apply Green's first identity 
and neglect the boundary integral. The result is $g ~ 
$si + *S2, where $si = (0i277i7i2 - 0^2"! • n2) /(I - 713) 
and $S2 = [0222772 - 20^22"2(n2 ■ n2)] /[47r(I - 773)^]. 
The coefficients 0i2 = 0i/(2a;), 0^2 = 0i/2, 0222 = 
(0i/2 + 202)/a;^ and 

0222 = <^'i/8 are functions of 
and the new vectorial weighted density is 



Rn2 



1 — n; 

ni(r) = -V772(r)/(8^). 

Replacing 771 only in half of the split term 77i0i/i? has 
two important consequences. Firstly, it insures that the 
function 0222 remains finite in the limit x — >■ 0. Sec- 
ondly, it leads to a symmetry between 0i2 and 05'2 that 
is reflected in their equality in the limits a; — > and 
a; 00, where 0i2(O) = 0i2(O) = —-^ and 0i2(oo) = 
05'2(oo) = \/2 — 2. However, comparing 0i2 with 05^2 
for any x reveals their differences (Fig. [T]). When we re- 
place 0^2 with 012, ^si is unaffected in the bulk fluid. 
Above we calculated the weighted densities . and 
in the OD limit. These are now supplemented by 
,^OD(r) = ijS'{R-r)/{STT), andnOD(r) = rjS' {R-r)r / {Stt) , 
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FIG. 1: CoefRcient functions in the attractive part $s of the 
excess free energy density for SHS with different Baxter pa- 
rameters r. In $3 the coefficients are functions of a; = , 
which in the homogeneous bulk fluid is related to the packing 
fraction rj via x = '^'^^ non-divergent expression $s, 

Eq. ([2]), relies on replacing by 012 and (j)222 by 0222 in "I>s. 



from which we have ni^n2^ 
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this modification regularizes $si 
vanishes as a;~^lna; fo 



■'12 
^X2 
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0. Hence, 
Note, that because 
r X — ^ oo then even if 



and 



>12 



take the same values in this limit it is insuf- 



ficient to insure that $gi vanishes in the OD limit, where 
large values of x occur due to the delta- function in ri2^. 
Furthermore, larger powers of the vectorial weighted den- 
sities could be used to modify the terms 0i2?t-i?t-2/ (I — Jt-s) 
and (/>i2ni • n2/(l — n^) individually while maintaining for 
example the PY direct correlation function which under- 
lies 4>g. However, the symmetry of (j)i2 and motivates 
the simple remedy introduced here. The value of giving 
up Baxter's PY result for the direct correlation function 
is assessed in terms of the performance of the modified 
functional in comparison to simulations. 

The regularization of $S2 requires two steps. We re- 
place 0222 by 0222 and note that even for asymptotic 



values of the argument x, 0222(0) 
'^'222(0) = - j^; and 0222(00) = 
while 0:^0, (cx)) = -i(2 - \/2) = 



]_ 

12t 648t^ 

i(V2-l)^- 
0.146, these function 



while 
-0.138 



and 

differ (see Fig.[T|). Therefore, our substitution is justified 
by the fact that it becomes exact in the limit of low den- 
sity {x 0) and weak "stickiness" (r 00). Like $si, 
<i>S2 in the bulk fluid remains unaffected by the substitu- 
tion. 
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FIG. 2: Pair-correlation function g{r) of a SHS fluid with 
r = 0.13 and rj = 0.164 from MC simulations (solid line) 
compared to results from the PY-based FMT (dotted line) 
and the regularized FMT (dashed line). 



After replacing 0222 by 0222 j ^S2 is still divergent. This 
is a consequence of the factor 2 in front of the vectorial 
term, which comes from matching the x — >■ and t — > 00 
limits of 0222 and 0222 ■ In analogy to $hs where the term 



~ 3n2(n2 • n2) is regularized by replacing it with n2 
we now regularize 712 — 2n2(n2 • n2) by substituting n2n|. 
Introducing the modified weighted density ni = ni — (ni • 
n2)/«2 (which equals ni in the bulk fluid and vanishes in 
the OD limit) we can write down the regularized version 
of <I>g as 



and 



512^1712 



£)222"-2'^2 



1 - ng 47r(l - nsY 
!)222 are functions of x 



(2) 



Rn2 
1— "3 ■ 



where 

As a first test of this regularized free energy density $3 
we calculate the pair correlation function g{r) by mini- 
mizing the density functional n = Jid + / dr[$HS + ^"3 + 
p(r)(T4xt(r) — /J.)], where J^id is the functional for the ideal 
gas, fi is the chemical potential and Voxt is the external 
potential, chosen so that it represents a particle of the 
SHS fluid at position r = 0. The resulting density pro- 
flic p(r) gives g(r) = p{r)/ p\j, where pb is the bulk fluid 
density. In Fig. [2] we show the case of strong adhesion, 
T = 0.13, and compare the FMT results with MC simula- 
tions obtained using the algorithm described in [ij] . The 
regularized 3>g is superior to the PY-based $Si which is 
remarkable because $s has been optimized for the CD 
limit and not at all for g{r). Other, larger r, MC simula- 
tions for g{r) (see [l^) are equally well described by the 
two FMTs, which give very similar results for r > 0.2. 

Alternatively, an analytical expression for the direct 
correlation function c(r) can be obtained from the excess 
free energy J"ox — J dr($Hs + $s) through 



c(r) 



Sp{r)6p{r') 



1^ ^o0{r), 



driadnp 



(3) 

where a and /3 run through {1, 2, 3, 1, 2} and r = |r — r'|. 
The partial derivatives of I'hs + are evaluated in 
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FIG. 3: The percolation line as obtained from MC simulations 
[l^ compared to the PY result (dashed line) and the results 
based on the non-diverging FMT constructed here (solid line) . 



the bulk fluid and derivatives with respect to vectorial 
quantities ni and n2 are executed formally while the 
fact that ni = and n2 = in bulk ensures that the 
resulting expressions are scalars. The symmetric co- 
efficients ujap are convolutions of the weight functions 
that arc associated with the respective weighted densi- 
ties. Therefore coap = for r > 2R. For example, — 
(4i? + r) {2R — r)^ is the overlap volume of two spheres 
with radius R and ccntcr-to-center distance r. More- 
over, UJ23 = ttR{2R - r), a;i3 = -r/8 + R/2 - i?V(4r), 
UJ22 = 27ri?Vr, UJ12 ^ R/{ir) + RS{2R - r)/8, CJ22 = 
2TTR^/r - irr, and W12 = i?/(4r) - R5{2R - r)/8 are 
required in Eq. ([3|). 

In particular, the coefficient of the (5-function in c(r) 
" + ■ This result can be used to 



IS a = — 



(l-r,) = 



calculate the percolation line which, under the approxi- 
mations described in 17[, is ^ = ry. As long as r is not 



too large, both the classic PY result and the modified 
FMT describe the percolation threshold as determined 
in simulations, the latter being somewhat more reliable 
in the range of intermediate r (Fig. [3]). 

Finally, we consider a confinement scenario for which 
our theory is expected to be the most relevant. We com- 
pare the previous PY-based FMT and the present reg- 
ularized FMT to canonical MC simulations for a SHS 
fluid with T = 0.13 and r = 0.2 in a spherical cavity of 
radius^pav ~ 3.7i?. This system was considered previ- 
ously [18[ for a hard-sphere fluid (t — > 00). To compare 
the FMTs with the simulations we have to transform the 
former from the grand canonical ensemble to the canon- 
ical ensemble which can be achieved up to corrections of 
order 1/A^'^ by using the scheme of [18[ that is based on 
results from [19[ . Fig. U shows that the regularized FMT 
is clearly superior to the previous PY-based FMT, espe- 
cially in the center of the spherical cavity. Interestingly, 
while the previous FMT underestimates the density in 
the center when r is small, it overestimates the density 
in the center when applied to hard spheres; large r (see 
[ist). The crossover occurs around r = 0.5 where the 



FIG. 4: Density profiles for A'^ = 10 sticky hard spheres with 
Baxter parameter r = 0.13 and r — 0.2 in a spherical cavity 
of radius 7?cav = 3.7R. The FMT results of the previous PY- 
based FMT and the present regularized FMT were obtained in 
the grand canonical ensemble and transformed to the canoni- 
cal ensemble (with correction of order 1 /N'^ ) for a comparison 
with canonical MC simulations that were performed using the 
algorithm of [l4| . 



regularized and the PY-based FMT yield very similar re- 
sults. In the limit of large cavity radius i?cav, planar wall 
density profiles are recovered; comparison (not shown) 
with existing MC simulations [2^ shows that the PY- 
based and the regularized FMT describe the density pro- 
files equally well. In particular, the theories yield identi- 
cal contact values pc for the density at a hard wall. This 
results from $hs+^s being equal to <E'HS+'i's in bulk and 
hence the underlying pressure being the PY compressibil- 
ity result ppy for both theories. Given that both FMTs 
obey the contact theorem ~ Ppvy the contact val- 
ues of the density must be identical. This illustrates the 
importance of preserving the bulk fluid properties while 
removing divergences for highly peaked density distribu- 
tions in order to obtain a density functional theory that 
is robust and accurate for a broad spectrum of settings. 
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